Experimental observation of optical rotation 
generated in vacuum by a magnetic field 
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Abstract 

We report the experimental observation of a light polarization rotation in vacuum in the presence 
of a transverse magnetic field. Assuming that data distribution is Gaussian, the average measured 
rotation is (3.9 ± 0.5) 10" 12 rad/pass, at 5 T with 44000 passes through aim long magnet, with 
A = 1064 nm. The relevance of this result in terms of the existence of a light, neutral, spin-zero 
particle is discussed. 
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Introduction - The existence of quantum fluctuations opens the possibility of considering 
vacuum itself as a material medium with dielectric properties. One such predicted property 
is the fact that vacuum becomes birefringent in the presence of a magnetic field due to 
photon-photon interactions 0,0,0]. Furthermore, the Primakoff effect, where light, neutral, 
scalar or pseudoscalar particles are produced from a two-photon vertex, could also give rise 
to macroscopically observable optical phenomena. In this second case it is expected that, 
by applying an external magnetic field, vacuum becomes birefringent and dichroic ODD 
The result of the induced dichroism is that the polarization plane of linearly polarized light 
is rotated by an angle 

(an — a_i_) „ , , 

a = K " D sin2fl (1) 

after traveling a distance D in a direction perpendicular to the field. Here a\\ and a± are the 
absorption coefficients of the parallel and perpendicular components of the light polarization 
with respect to the magnetic field, and 9 is the angle between the magnetic field and the 
polarization vector. 

The PVLAS collaboration is studying the magnetic birefringence and dichroism of vac- 
uum with a high-sensitivity (~ 1CT 7 rad/yTJz) optical ellipsometer, at the INFN Legnaro 
National Laboratory, Legnaro, Italy 0. In this apparatus, a linearly polarized light beam 
at 1064 nm from a Nd:YAG laser interacts, inside a Fabry-Perot optical resonator, with the 
field produced by a superconducting, 1 m long dipole magnet operated at a maximum field 
of 5.5 T. 

In this letter we report the first observation of a rotation of the polarization plane induced 
in vacuum on linearly polarized light by a transverse magnetic field. This signal is generated 
within the Fabry-Perot cavity when the magnet is energized. Data from these measurements 
can be interpreted in terms of the existence of a very light, neutral, spin-zero particle coupled 
to two photons. 

Apparatus and experimental technique - Figure 1 shows a schematic layout of the PVLAS 
apparatus. The 1 m long interaction region is contained within a 4.6 m long, 25 mm diameter 
quartz tube. This tube is placed vertically and traverses the bore of the dipole magnet. The 
magnet is housed within a 5 ton, 3.1 m high, warm-bore cryostat, and the field direction 
lies in the horizontal plane: the cryostat itself sits on a 0.9 m radius turntable standing on 
a concrete beam fastened to the experimental hall floor [8j. The UHV vacuum chambers 
hosting the optical elements, placed at the top and at the bottom of the tube, are fastened 
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Figure 1: Schematic layout of the PVLAS experimental apparatus. PI and P2 crossed polar- 
izing prisms, Ml and M2 Fabry-Perot cavity mirrors, QWP quarter-wave plate, SOM ellipticity 
modulator, RGA residual gas analyzer. 

to two granite optical benches. These benches are fixed on a concrete platform mechanically 
isolated from the rest of the hall floor. During normal operation the cryostat is filled with 
liquid He at 4.2 K and the magnet is energized with a current of about 2030 A, resulting in 
a 5.5 T field over the entire interaction region. To allow rotation of the magnet, the coils are 
shorted and disconnected from the power supply. The field intensity is monitored by a set of 
Hall probes [8]. The turntable is actuated by a low- vibration hydraulic drive and normally 
rotates the magnet-cryostat assembly, around a vertical axis, at a frequency v m ~ 0.3 Hz. 
During vacuum measurements, the quartz tube and the two vacuum chambers are kept 
in UHV conditions (P ~ 10~ 8 mbar) by a liquid N 2 trap combined with Ti sublimation 
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pumps. This pumping scheme has been chosen in order to avoid mechanical vibrations and 
possible couplings between the rotating dipole field and the ion-pump permanent magnets. 
The residual gas composition is monitored by means of a residual gas analyzer (RGA). 
The interaction region is contained within a high-finesse Fabry-Perot optical resonator (FP) 
formed by a pair of dielectric, multilayer, high- reflectivity, 11 m curvature radius mirrors 
placed 6.4 m apart (Ml and M2). Resonator diffraction losses are found to be negligible. 
The ellipsometer consists of a pair of crossed polarizing prisms (PI and P2) together with 
an ellipticity modulator (Stress Optic Modulator, or SOM The SOM, driven by a Sine 
Generator at a frequency vsom = 506 Hz, provides a carrier ellipticity which beats with the 
modulated signal produced by magnet rotation and allows heterodyne detection. A quarter- 
wave plate (QWP), inserted between the M2 mirror and the SOM, can be placed in or 
removed from the beam path without breaking the vacuum. It can also be rotated in order 
to exchange its fast and slow axes. The laser beam (from a 100 mW, CW, Nd:YAG laser at 
.064 nm) is frequency-locked to the FP cavity by means of an electro-optical feedback loop 
. In this way, the optical path is increased by a factor N = 2JF/7T, where T is the finesse 
of the FP. Typically, for the PVLAS cavity, T ~ 7 10 4 . Light transmitted through the 
analyzer P2 is detected by a photodiode followed by a high gain (10 7 V/A) transimpedance 
amplifier, and the resulting voltage signal is fed to the digitizing electronics. The position 
of the rotating table is sensed by an optoelectronic position detector and this signal is also 
read and digitized: this allows determination of absolute phases. 

Measurements and results - The light intensity transmitted through the crossed polarizers 
of the PVLAS ellipsometer is [3] 

/ = /oK + Kt) + rKt) + r(t)] 2 } (2) 

where Iq is the intensity before the analyzer P2, a 2 is the extinction factor, 
a(t) = a cos(47rz/ m £ + 20 m ) is the rotation to be detected, r](t) = r] Q cos(27tz/ S om^ + #som) is 
the carrier ellipticity, and T(t) represents additional, quasi-static, uncompensated rotations 
and ellipticities. Table I lists the Fourier components obtained by expanding expression |2l 
and it shows that the amplitude of the two sidebands separated by 2z/ m (twice the magnet 
rotation frequency) from the carrier frequency at z^som is proportional to the the amplitude 
of the rotation to be detected. When properly aligned, the QWP transforms apparent ro- 
tations (dichroisms) into ellipticities [lOj, which then beat with the carrier ellipticity signal 
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Frequency Intensity /Iq Phase 

a 2 + T 2 + r? 2 /2+a2/2 
4^ m a 2 /2 40 m 

^SOM 2r?7o #SOM 

z-'SOM ± 2z/ m 770^0 #som ± 2# m 

2^som ^o/2 26>som 

Table I: Main frequency components of the detection photodiode signal. 

and are detected. These ellipticities conserve amplitude and undergo a sign change if the 
fast axis of the QWP is at 90° with respect to the initial polarization direction. We label 
this setting as QWP 90°, and use QWP 0° when the fast axis of the QWP is aligned with 
the initial polarization direction. In both cases, ellipticities are converted into rotations and 
pass undetected. Conversely, without the QWP in the beam path, ellipticities acquired in 
the interaction region beat with the SOM carrier ellipticity and are detected, while acquired 
rotations do not beat with the SOM, and pass undetected. 

As a check of the Fourier phase of observed physical signals, we have measured the Cotton- 
Mouton effect (CME) [12] in N 2 and Ne at several pressure values (for N 2 pressure ranged 
from 6 /ibar to 1.8 mbar, and for Ne from 0.5 mbar to 20 mbar). For these measurements, 
the gas was inserted in the interaction region, and the QWP was momentarily removed. 
Sideband peaks were then observed in the photodiode signal Fourier spectrum at twice the 
magnet rotation frequency from the carrier. The phases of these peaks correspond to the 
angular position of the turntable at which the field direction is at 45° with respect to the 
initial fixed polarization. In this position the magnetic birefringence due to the CME is 
maximum, and the numerical value of the phase is fixed by the geometry of the apparatus. 
Phase and amplitude data resulting from these measurements were plotted on a polar plane 
and fitted with radial lines. These lines fix what we call the physical axis, and are in 
agreement with geometrical considerations. 

Two typical Fourier amplitude spectra of the detection photodiode signal in vacuum are 
shown in Figure 2. Figure 2a shows a spectrum around the 506 Hz carrier frequency when the 
magnet is off: this spectrum does not show sideband peaks at ±2i/ m from the carrier. Figure 
2b shows an amplitude spectrum observed in the same conditions of Figure 2a, but with the 
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Figure 2: Typical Fourier amplitude spectra from vacuum rotation data (P ~ 10 -8 mbar) with and 
without magnetic field. Both spectra were taken with the magnet rotating. Arrows and numbers 
below the curves indicate the expected position of sidebands with frequency shifts that are integer 
multiples of u m ; the relevant signal peaks correspond to frequency shifts of ±2u m . 



magnet on at 5.25 T. Here sideband peaks at ±2z/ m are present (peak phase is 336° ± 6°). 
To simplify the comparison, we show spectra with a similar noise level (~ 10~ 8 rad) near the 
"2" sidebands. The sideband peaks at frequency shifts ±i/ m were studied and found to be 
uncorrelated, both in amplitude and in phase, to the signal sidebands at ±2z/ m . Sideband 
peaks at a distance ±2z/ m from the carrier, similar to the ones appearing in Figure 2b, 
have been observed in all the PVLAS data runs with the magnet energized, and sideband 
peak phases have been measured in addition to their amplitudes. Several tests have been 
performed to clarify the nature of these peaks. We find that the peaks are present only when 
the magnet is on and the mirrors of the FP cavity are mounted. We have also observed 
that when QWP is removed from the beam the peak amplitude changes. Further, when the 
fast and the slow optical axes of the QWP are exchanged, that is the QWP is rotated by 
90°, signal phase changes by 180° as expected for a "true" rotation Q], while the amplitude 
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Figure 3: Polar plot showing the distribution of rotation data in vacuum with B = 5 T and (a) the 
QWP in the 90° position, and (b) the QWP in the 0° position. Each data point corresponds to a 
single 100 s record. The statistical uncertainties from the fit procedure are small (~ 10~ 12 rad/pass) 
and are omitted for clarity. Solid lines represent the resulting average vectors and the dotted line 
the physical axis. 

remains unchanged. These facts indicate that the "±2i/ m " sidebands are associated with an 
interaction inside the FP cavity. 

Vacuum data have been collected for two years, in the course of six two-week runs, and 
were routinely taken with cavity finesses greater than 62000. Collected data were subdivided 
in 100 s long records, for a total of 84 records at QWP 90° and 121 records at QWP 0°. Each 
record was then analyzed with a fitting procedure that compensates the mechanical jitter 
of the rotating table and the uneven sampling of the table position Q|. For each record, 
the procedure returns the x and y components of a vector in the amplitude-phase plane of 
the "±2z/ m " sidebands; these components are Gaussian variates, and the procedure returns 
their statistical variances as well. Since these data records have been taken with different 
cavity finesses, to plot all available rotation data records in a single graph amplitudes have 
been normalized to a single pass in the cavity. The measurements performed with the two 
different QWP orientations, taken within the same two-week run, appear in the x — y plane 
as two groups of points located on opposite sides with respect to the origin. However, the 
average phase and amplitude of each group vary from one run to the other more than their 
individual dispersion. As a consequence, the standard deviation of the entire set of x and y 
components shown in Figure 3 is larger than the standard deviation obtained from the fit 
procedure: the values are (a) ((Js,x, &s,y) = (2.3,1.8) 10~ 12 rad/pass, and (b) (as,x, vs,y) = 
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Table II: Values and corresponding standard deviations, both statistical (stat) and total (tot), in 
units of 10 -12 rad/pass, of the x and y components of the vectors shown in Fig. 4. 

(2.4,2.2) 1CT 12 rad/pass. We interpret this as an indication of the presence of unidentified 
systematic effects which increase the overall measurement uncertainty. Assuming that these 
systematic effects have a Gaussian distribution in the long term, we estimate additional 
variances o~g x , a 2 Sy from the distributions of the x and y components. Thus we associate 
to each x and y component a new variance o\ ot which is the sum of the statistical variance 
from the fit procedure and of the additional variances a\ and a| , and calculate weighted 
averages of x and y for each orientation of the QWP. Figure 4 shows both the weighted 
average vectors, and 90, and the half-difference vector A, representing that part of the 
signal changing sign under a QWP axis exchange. The half-sum vector £ is a systematic 
error. This analysis suggests a physical origin for the rotation signal observed in vacuum 
when the magnetic field is present (numerical values for X and A are listed in Table 11). 
From Table II we find the following value for the amplitude of the measured rotation in 
vacuum with B m 5 T (quoted with a 3<r uncertainty interval): 

a = (3.9 ± 0.5) 10~ 12 rad/pass (3) 

Discussion - We have observed a rotation of the polarization plane of light propagating 
through a transverse magnetic field having the following characteristics: 
-it depends on the presence of the magnetic field; 
- it is localized within the Fabry-Perot cavity; 

-it has the proper phase with respect to the magnetic field instantaneous direction. 

Experimentally, then, we find that vacuum induces a rotation of the polarization plane in 
the presence of a magnetic field. The possibility that this effect is due to an unknown, albeit 
very subtle, instrumental artifact has been investigated at length without success. A brief 
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Figure 4: Polar plot showing the weighted averages of the vectors shown in Figure 3. Angles are 
measured in degrees and amplitudes in rad/pass. The dotted line at 15.1° shows the direction of 
the physical axis. The vectors and 90 represent the average rotations for two orientations of the 
QWP, with ellipses at their tips giving the 3a uncertainty regions. The vectors A and £ are the 
half-difference and half-sum of the average vectors, respectively. 

listing of the possible spurious effects which have been considered includes: electrical pick- 
ups, effects from the residual gas, beam movements induced by couplings of the rotating field 
to parts of the apparatus (the permanents magnets of the ion vacuum pumps, for instance), 
higher harmonic generation by the sideband peaks at ±z/ m , direct effect of the field on the 
optics. 

Speculations on the physical origin of the measured signal suggest the possible presence 

of new physics. In fact, a selective absorption of photons, resulting in an apparent rotation, 

n 

could be caused by photon splitting processes |2[, where photons interact with the external 
field and split into two or more lower energy photons. These photons would not resonate with 
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the Fabry-Perot cavity and would go undetected through the PVLAS ellipsometer, where 
this effect would be seen as an absorption. However, in our apparatus photon splitting is 
shown to be negligible. 

Another source of apparent rotation could be associated to a neutral, light boson produced 
by a two-photon vertex P, 0, Q] . In this photon-boson oscillation framework, both induced 
rotation and induced ellipticity can be expressed as functions of the particle mass m& and 
of the inverse coupling constant to two photons Mf, . In particular, the amplitude of the 
induced dichroism e is given by (in natural Heaviside-Lorentz units) 



sin 29 } X 

\AM b 



r sin (m\L/A.uj) 



(4) 



where 9 is the angle between the magnetic field B and the light polarization direction, 
L is the length of the magnetic region, N is the number of passes in the FP cavity, and 
uj is the photon energy. Therefore, the measured rotation in vacuum fixes a relationship 



between mass nib and inverse coupling Mf,. It is important to note that if we set in formula 

-i 2 

s'm(m?L/4uj) — . . 

= 1, the constant Mf, implied by result |3| is bound 



0] the diffraction term 



by Mi % 6 10 5 GeV. This is compatible with rotation limits published in When 
these limits are combined with the Mb vs. nib curve corresponding to the PVLAS signal, we 
find that only one segment of this curve is allowed as a region for possible (nib, Mb) pairs. 
The allowed segment is contained in the intervals 2 10 5 GeV < M b ^ 6 10 5 GeV, and 
1 meV < m 6 < 1.5 meV. 

Measurements are already under way with a slightly modified version of the apparatus, 
where a Faraday-cell type polarization modulator has been inserted in the beam path to 
obtain simultaneous detection of rotations and ellipticities. During these measurements, 
Ne gas has been inserted at several pressures in the interaction region in order to change 
the optical path: this will allow a direct test of the supposed photon-boson oscillation 



phenomenon 



12|. 



Recently, limits on axion-photon coupling with a strength far from this estimate have 
been published in While a detailed discussion of astrophysical bounds is beyond the 
scope of the present letter, we mention that recent theoretical work (l^| suggests scenarios 
where it is possible to accommodate both results. 

Let us finally note here, that a crucial check of these speculations could come from a 
"photon- regeneration" experiment of the type already attempted in and in jsoj l . 

10 



Acknowledgements - We wish to thank S. Marigo and A. Zanetti for their invaluable 
and precious technical help on the construction and running of the apparatus, and also 
R. Pengo and G. Petrucci who designed and built the rotating magnet system. We also 
acknowledge all the people of the Laboratori Nazionali di Legnaro that helped us and made 
this experiment possible. 



[1] W. Heisenberg and H. Euler, Z. Phys. 98, 718 (1936); V.S. Weisskopf, Mat.-Phys. Medd.-K. 

Dan. Vidensk. Selsk. 14, 6 (1936); J. Schwinger, Phys. Rev. 82, 664 (1951). 
[2] S.L. Adler, Ann. Phys. 67, 559 (1971). 
[3] E. Iacopini and E. Zavattini, Phys. Lett. B 85, 151 (1979). 

[4] P. Sikivie, Phys. Rev. Lett. 51, 1415 (1983); A. A. Anselm, Yad. Fiz. 42, 1480 (1985); M. 

Gasperini, Phys. Rev. Lett. 59, 396 (1987). 
[5] L. Maiani, R. Petronzio and E. Zavattini, Phys. Lett. B 175, 359 (1986). 
[6] G. Raffelt and L. Stodolsky, Phys. Rev D 37, 1237 (1988). 

[7] D. Bakalov et al, Nucl. Phys. B (Proc. Suppl.) 35, 180 (1994); D. Bakalov et al, Quantum 

Semiclass. Opt. 10, 239 (1998); D. Bakalov et al, Hyperfine Int. 114, 103 (1998). 
[8] R. Pengo et al, in Frontier tests of QED and Physics of the Vacuum, E. Zavattini et al. eds., 

Heron Press, Sofia, 59 (1998). 
[9] F. Brandi, E. Polacco and G. Ruoso, Meas. Sci. Technol. 12, 1503 (2001). 
[10] See for instance M. Born and E. Wolf, Principles of Optics, Pergamon Press, 691-692 (1980). 
[11] G. Cantatore et al., Rev. Sci. Instrum. 66, 2785 (1995). 

[12] S. Carusotto et al, J. Opt. Soc. Am. B 1, 635 (1984); C. Rizzo, A. Rizzo and D. M. Bishop, 
Int. Rev. Phys. Chem. 16, 81 (1997); M. Bregant et al, Chem. Phys. Lett. 410, 288 (2005). 
[13] E. Milotti, J. Comput. Phys. 202, 134 (2005). 
[14] E. Masso and R. Toldra, Phys. Rev. D 52, 1755 (1995). 

[15] In this model, the induced rotation is caused by real production of bosons, while the induced 

ellipticity is due to virtual production. 
[16] R. Cameron et al, Phys. Rev. D 47, 3707 (1993). 

[17] E. Zavattini et al, in 11th International Workshop on Neutrino Telescopes, M. Baldo Ceolin 
ed., Papergraf- Padova (Italy), 433 (2005). 



11 



[18] K. Zioutas et al, Phys. Rev. Lett. 94, 121301 (2005) 
[19] E. Masso and J. Redondo, |a7Xiv:hep-ph/0504202| 
[20] G. Ruoso et al, Zeit. f. Phys. C 56, 505 (1992). 



12 



